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ISOSINGLET APPROXIMATION FOR NONELASTIC REACTIONS 


By John W. Wilson 
Langley Research Center 

SUMMARY 

It is often necessary in high-energy transport calculations to approximate in some 
manner the cross sections for which data are not available. Most common are interpo- 
lation in energy or atomic number and mass. Relations of yet another kind can be estab- 
lished by using isospin symmetries of the strong force and a few simple assumptions 
about nuclear reactions and nuclear states. 

The Serber transition operator is formed of two factors t = t g t c , the cascade and 
evaporation part, respectively. The cascade event is via strong interactions which con- 
serve isospin, and the evaporation is isospin symmetry breaking. 

Since t c proceeds via the strong force, then isospin rotation relates cascade 
cross sections of isospin conjugate states. That isospin conjugate states of the evapora- 
tion cross sections are not similarly related follows since T e is isospin symmetry 
breaking. Since the strong force is most cohesive for states of lowest isospin, however, 
it is reasonable to assume that the intermediate nuclear state R is an isosinglet (that 
is, an isospin scalar particle). Thus, the state R is unchanged under isospin rotation 
and the evaporation spectra are consequently unaffected by isospin rotations. 

INTRODUCTION 

High-energy hadron transport calculations require an extensive knowledge of the 
nonelastic cross sections for processes occurring within the transport media. Although 
considerable progress has been made toward attaining complete sets of experimental and 
theoretical data, they are as yet very incomplete. In order that transport calculations 
could proceed, approximate relations have often been developed to extend existing data 
into regions where cross sections are unknown (refs. 1, 2, 3, and 4). 

It is well known that energetic nuclear production cross sections are smooth func- 
tions of A and Z for sufficiently large A. Standard procedures have developed which 
use this property in preparation of cross sections for transport calculations (ref. 1). This 
procedure has eliminated the need for performing experiments for many materials for the 
purpose of nucleon transport studies. 



Methods have been developed for interpolation and extrapolation in primary energy 
(refs. 1, 3, and 4). They are based on the observation that the high-energy secondary 
production spectra, if plotted against the ratio of secondary energy to primary energy, 
are found to be nearly independent of primary energy (ref. 1). Although this property is 
expected to be true only between inelastic thresholds (refs. 2 and 3), detailed studies have 
shown it to be not unreasonable in extrapolating to very high energies over many such 
thresholds (ref. 4). Caution must be applied with this method since it applies only to the 
so-called "cascade" event and is not valid : for the low-energy nucleons usually attributed 
to "evaporation" processes (refs. 1 and 5). 

The present work will describe in detail extensions of yet another type. It gener- 
ally occurs that experimental data are available for limited types of primary and secon- 
dary particles; for example, high-energy beams of neutrons are rarely available, and 
neutrons are not usually measured reliably in the final state. Thus' neutron data almost 
always must be inferred from measurements on other particles. Described herein is a 
derivation of relations based on an SU(2) approximation of the nonelastic process; 

These relations have proven to be helpful in developing cross -section data for transport 
calculations (refs. 2, 3, and 5). The development uses the Serber two-step model (ref; 6), 
but the results may have validity in a broader sense. 

SYMBOLS 

A nucleon number 

E entering-state label ... 

F final-state label ■ v • 

F’ transformed final -state label , ti ,• ... . i .„, 

F c final cascade -particle -state label 

F e , final evaporation-state label . f . ...... 

Fj intermediate scattering -state label 

ly(0) " isospin operator for rotation through angle 9 about iy 

i,i z ,i’,i z isospin quantum numbers, dimensionless 
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iy unit vector along y-axis in isospin space 

N neutron number 

P proton number 

p,n nucleon-state labels, p for proton state and n for neutron state 

R residual excited -nuclear -state label 

T transition-operator matrix elements, amu _ l 

u general isospin rotation operator, dimensionless 

K,fi energy and angle of secondary particles 

6jj Kronecker 6 -function 

r / phase factor 

n^n 0 pion-state labels (superscript denotes charge) 

t state transition operator, amu'^ 

| a) vector in Hilbert space with label a that represents a complete set of 

quantum numbers 

CT X Y double differential nonelastic cross section from entering state X to final 
state Y 

or Y v double differential nonelastic cross section for production of y -particle with 

A ? y 

x -projectile 

Subscripts and superscripts: 
c cascade 

E entering 
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e 


evaporation 


F final 

Fj intermediate state 

p projectile 

R residual 

s strong 

y y -component 

z z -component 

y electromagnetic 

t dagger over an operator denotes Hermitian conjugate 

ELEMENTARY CONCEPTS 
Serber Model 

It will be assumed for the present purposes that the Serber two-step model (ref. 6) 
for high-energy nuclear reactions is substantially correct. The dynamical basis of this 
model has been discussed extensively elsewhere (refs. 6, 7, and 8 and references therein) 
and will not be repeated here. The cascade process of the first step proceeds via the 
strong interaction (the coulomb scattering peak is neglected) and is subject to the sym- 
metries (in particular, isospin) of the strong force. The first-step transition operator r c 
is defined so that the transition matrix element is 

T c (Fi,E)= <Fi|t c |E) (1) 

in state vector notation (ref. 9) where |E) and |Fj) denote the asymptotic entering 
and final states of the cascade process, respectively. In the model, it is assumed that 
the final cascade state contains one or more highly excited residual nuclei denoted by 

l F l> = l R >l F c> 
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where |R) is the residual nuclear state and |F C ) is the cascade particle state (that is, 
particles ejected by direct interaction). The excited nuclei decay through particle emis- 
sion (evaporation) and finally y -radiation. This deexcitation occurs through strong and 
electromagnetic interaction which breaks isospin invariance. The step -two deexcitation 
transition operator r e is related to the total T -matrix element of the reaction by 

T= <F|T e rc|E>=^<P|r e |F 1 >(F, 1 |r c |E> (2) 

Since the Serber model conveniently divides into these two steps, one which is iso- 
spin invariant and the other which is isospin symmetry breaking, this reaction is easily 
analyzed. 


Isospin, SU(2) 

An outline of those properties of the representations of the isospin group in rela- 
tion to observations of hadronic systems which are required is now presented. Addi- 
tional detail on SU(2) and the strong interaction are to be found in reference 10. Let 
u e SU(2); then, for every transition operator r s of only strongly interaction particles 



which is the statement of isospin invariance. The quantum numbers of the isospin group, 
i and i z , are conserved by r s as follows: 

<i',i z |T s |i,i z > = T s 6 ii ,6 i z iz ( 4) 

where superscript i on T s means it depends only on i and not on i z . Electromag- 
netic interactions r y depend on the projection i z , and it follows that generally 



and isospin is not conserved. 

Systems of hadrons are identified by isospin quantum numbers which belong 
to various isospin multiplets. The multiplicity determines the i state (that is, 
Multiplicity = 2i + 1 for i = 0, 1/2, 1, 3/2, . . .) with the i z component measured 
by the symmetry -breaking electromagnetic interaction (that is, electric charge). The 
symmetry operation of importance here is isospin rotations about the iy axis by n 
as denoted by Iy(7r). The following hadrons are of interest: 
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(1) Nucleon isodoublet 

There are two change states: 


Proton |p) = 1 1/2, 1/2) (6a) 


Neutron |n) = 1 1/2, -1/2) 

(6b) 

which are related by 


£ 

ii 

(6c) 

Iy(77)|n) = - |p) 

(6d) 

(2) Pion isotriplet 


There are three pion-charge states: 


k + > = Im> 

(7a) 

=3 

^o. 

II 

(7b) 

k"> 

(7c) 

which satisfy 


lyOOU*) = -U T > 

(7d) 

£ 

ii 

i 

=i 

(7e) 

(3) Nuclei 


Nuclei which are presently considered to be bound systems of nucleons are char 
acterized by nucleon number A, proton number P, and neutron number N where 

A = P + N 

(8a) 

i z = (P - N)/2 

(8b) 

It is observed that small isospin numbers are favored by the 
ground-state nuclei have the smallest i consistent with i z 

strong force; for example, 
which gives 
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i = |i z ) = |P - N|/2 


(8c) 


The operation Iy(tf) yields 

IyOOlMz) = *?|Wz) ( 8d ) 

where 77 is a trivial phase factor. For nuclei with P = A/2 

i = |P - N|/2 = |A - A |/4 = 0 (8e) 

which are isosinglet hadrons. 

■•o 

ISOSINGLET APPROXIMATION 
Step One, Cascade 

The cascade process consists of multiples of binary collisions within the target 
nucleus. The number of binary collisions with nucleons of the same type clearly depends 
on the number of nucleons of that type (that is, P or N, whichever is appropriate). The 
fractional number of collisions with protons is expected to be ~P/A and similarly for 
neutrons «N/A. Furthermore, if P and N are large, then (P ± 1)/A ~ P/A and 
(N ± 1)/A * N/A so that secondary production spectra are not expected to be rapidly 
varying functions of P and N. On this basis, the amplitudes for 

p + Anything 

p + C 13 — Anything 

differ by no more than 

(7/13 - 6/12) ~ 8 percent 

for cascade processes. In the following it is assumed that the neighboring nuclei com- 
posed of combinations of P ± 1 and N ± 1 are indistinguishable in nonelastic reactions. 
Furthermore, only nuclei for which 


P fc N « A/2 (9) 

are considered and refer to the isosinglet approximation as the assumption that 

P = N = A/2 (10) 
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The entering state consists of two hadrons, one of which has A » 1 for which the 
isosinglet approximation is assumed to be valid; so the isospin quantum numbers of the 
entering state are those of the second hadron (hereafter referred to as the projectile and 
denoted by subscripts p). Hence, using the isosinglet approximation, for the entering 
state 


IWE^E) = ||P-N|/2,(P-N)/2>|ip,ip z > 

~ l°»°)l ip^Pz;) = | i p» i p z , E > (H) 

where isospin numbers of the states are now written explicitly. 

The final state consists of the many hadrons which were expelled from the nucleus 
by the direct binary reactions of the cascade process and one or more heavy residual 
nuclei which are left in highly excited states. These excited nuclei subsequently decay 
through particle and y emission. These residual nuclei, it is assumed have 

p R « % * Ar/ 2 (12a) 

A r » 1 (12b) 

and, since the strong force favors states of small isospin, for residual nuclei 

iR=l p R-N R l/2”° (12c) 

that is, an isosinglet approximation. The cascade particles (those expelled by binary 
collisions) have total isospin number i c and projection i c . Hence, the final cascade- 
state configuration is 

l iF T iF l z ,F l) = l i R’ i R z ’ R )l ic ’ ic z ,Fc ) 

*|0,0,R>|i c ,ic z ,F c ) Hicic^Fj) ( 13 ) 

The assumed isospin invariance of the cascade transition operator implies for 
U = Iy(7T) 

T C ( F 1,E) = ( F ljicAc z | T c| *p>*p z >E) 

= < F l> i c. i c z l u+UT c utu i i p> i p z ,E) 

= ( F l> i c>~ i c z | T c| i p>' i p z >E) = T c P 6i pi c 5jp i c (14) 

^ z z 
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since r c = ut c u + , and in the last line of equation (14) u is taken as I y (7r). Hence, the 
cascade transition matrix element depends only on total isospin of the projectile, and the 
final-cascade -particles isospin is equal to the projectile isospin. 

Step Two, Evaporation 

As noted earlier, the evaporation process is greatly dependent on electromagnetic 
forces which are not isospin invariant. The evaporation characteristic time is long 
(»10 - ^3 sec ) compared with the short time duration of the cascade event (=10~23 sec); 
hence, the excited residual nuclear states are comparatively stable and the residual 
nucleus at the end of the cascade |R) is far removed in time from the cascade inter- 
action region. Thus, nontrivial matrix elements of the evaporation operator are taken 
between the excited residual nuclear state | R) and the final evaporation products | F e )> 
as 


< F e! T el R > 

The final state is naturally divided into two sets of particles: the cascade particles pro- 
duced in direct binary reactions and the evaporation products of the residual excited nuclei 
which are denoted as follows 

|F e > |K C > = |F> 

The full transition matrix element is 

T(F,E) = <F,i F ,i Fz |T e T c |i p ,ip z ,E> 

= J <F,iF>iF z | T e |i> i z , F l)( F l>Mz| r c | i p ,ip z> E) 

R i, Mz 

= £ <F,iF,iF z |T e | i p, i p z , F 1 ) T c P ( F l» E ) 

F 1 


= ^< F eKe|R)Tc P ( F l,E) (15) 

R 

where | F g ) is the state of final evaporation products and j R) is the state of residual 
nuclei appearing in the final cascade state and is assumed to be an isosinglet. The 
assumption that | R) is an isosinglet effectively decouples the element (F e jT e |R) in 
isospin rotations. 
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FINAL RESULTS 


It has been suggested that the nonelastic reaction is independent of the isospin of 
the target nucleus when 

!p-n|/a«i (16) 

so the entering-state isospin quantum numbers can be taken as those of the projectile. 
Similarly it is argued that the residual nuclei after the cascade processes can be approx- 
imated by isosinglet states. The final state is thus divided into two sets of particles: the 
cascade particles produced in direct binary reactions and the evaporation products of the 
residual excited nuclei which are denoted as follows 

|F e >l F c ) = |F) (17) 

Assuming the residual nuclei to be isosinglet effectively decouples the evaporation pro- 
ducts from isospin rotations u as seen by combining equations (14) and (15) as follows 

T(F\uE) =^(F e |r e |R)Tc P (uF 1 ,uE) (18) 

R 

where |F') = j F e )| uF c ) which follows from the assumption that |R) is an isosinglet, 
so that 

u|R> = |R) (19) 

Again the evaporation process due to assumption (19) has been effectively decoupled in 
isospin rotations. 

The cross section of these nonelastic processes is proportional to the T-matrix 
element squared. Equation (18) implies that the matrix element for the process E — F 
is numerically equal to the matrix element for (uE) — F'. The assumed invariance of 
the cascade operator and the isosinglet approximation imply that 

°E— F = ct (uE)— F* 

where |F’)= |F e )| uF c ). These results are applied to specific scattering systems in the 
next section. 
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APPLICATION 


Isosinglet -Isosinglet Reactions 
When both hadrons of the entering state are isosinglets 

u|E> = |uE> = jE) (20) 

ic ~ip = 0 (21) 

Thus, from equation (18) (F' as in eq. (18)) 

T(F*,uE) - V <F e |T e jR)T°(uF 1 ,E) (22) 

R 

from which it may be deduced that the cascade particles which are related by 

lyWlMz) = v\h~^z) (23) 


have the same energy spectra. For reactions induced by deuterons or a -particles the 
following correspondence between measured parameters of secondary cascade particles 
may be noted: 

cC^K.SJ) = 

^, 3 H (K ' n) = 

These relations are especially important for high-energy secondary neutrons since they 
are rarely measured reliably. 

Isosinglet -Isodoublet Reactions 

Among the isodoublet hadrons which are used as projectiles are (p,n) and (^He,^H) 
corresponding to i z = (1/2, -1/2). The two possible entering states are denoted by lE + ) 
and | E“) in an obvious notation with 

I y (7T)|E-> = p|E + > 
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where ri is an unimportant phase factor assumed to be 1. Applying equation (18) (for 

u = Iy(7T)) 


T(F’,E") =y <F e jT e |R)Tj /2 (uF 1 ,uE + ) 

R 

with F' as before. Hence the relations among secondary cascade particles for the two 
entering isospin states follow: 


''p, P ( K > n > = ^, n < K . n ) 
<£„<&>«> - o£ ;P <K,!2) 

a= T± (K,n) = d £ ;lrT <K,Si) 

a c 3 „(K,n) = a c 3u (K,SJ) 

p, J H n, J He 

ag >cl (K.n) = aS jd (K,iJ) 


pJ,a<K,n) = <,«(K ,n) 

which are production cross sections for protons, neutrons, deuterons, tritons, helion, 
a -particles, and pions due to incident protons and neutrons. Since the intermediate 
residual excited nuclear states are assumed to be isosinglet nuclei, the following rela- 
tions among secondary evaporation products are also, obtained: 


4 ,p<K.«> - <,p< K . n > 

p|„(K,SI) = c-e„(K,n) 
a* d (K,!2) = o« d (K,a) 


Similar results can be obtained for various entering -state configurations. 



Comparison of Results 


There are as yet insufficient experimental results to make comparisons with the 
present work. The only comparisons to be made are relations among results of model 
calculations of nonelastic reactions. Currently, the most complete model of nonelastic 
reactions is that of Bertini (ref. 11). The essential feature of the model of interest here 
is that it is a Serber model with isospin invariant reactions in the first step followed by 
a standard evaporation model (Dresner, ref. 12) for the second step. Extensive calcula- 
tions have been made for reactions induced by p, n, ir ± , and oversimplified a -particles. 
The cascade secondary particles allowed in the Bertini calculations (ref. 11) are p, n, 

■n ± , and 7T°. According to the present work, the high-energy proton and neutron spectra 
for a -induced reactions should be equal in the isosinglet approximation. Results obtained 
from the o: model of reference 13 (Gabriel and others) for reactions in ^Al using the 
Bertini calculations (ref. 11) are shown in figure 1. These results agree to within the 
statistical accuracy of the Monte Carlo method. 



E, MeV 

27 

Figure 1.- Comparison of secondary-production spectra for the a- 'A1 model reaction of 
reference 13. These spectra should be equal under the isosinglet approximation. 
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TABLE I. - QUANTITIES WHICH ARE EXACTLY RELATED BY SU(2) IN 
NONELASTIC REACTION. (DATA TAKEN FROM REF. 11.) 


Cascade multiplicity 

P_16q reaction 

n -16o reaction 

Proton 2.051 

Neutron 1.471 

?r+ 0.268 

v- 0.081 

77° 0.180 

Neutron 2.045 

Proton 1.490 

7 r- 0.280 

7 7 + 0.080,, 

7?o 0.165 


TABLE II.- QUANTITIES RELATED THROUGH THE ISOSINGLET APPROXIMATION 
APPLIED TO ENTERING STATE. (DATA TAKEN FROM REF. 11.) 


Cascade multiplicity 

p-^Al reaction 

n-^Al reaction 

Proton 2.291 

Neutron 1.861 

77+ 0.263 

77 - 0.087 

77 ° . 0. 165 

Neutron 2.426 

Proton 1.695 

77 - 0.273 

77 + 0.080 

77 ° 0.160 


TABLE III. - QUANTITIES RELATED THROUGH THE ISOSINGLET APPROXIMATION 
APPLIED TO THE INTERMEDIATE RESIDUAL EXCITED NUCLEAR STATE. 

(DATA TAKEN FROM REF. 11.) 


Evaporation multiplicity 

P-16q reaction 

n-l^o reaction 

Proton 0.841 

Neutron 1.114 

Deuteron 0.091 

Triton 0.025 

Helion ( 3 He) 0.029 

Alpha 0.341 

Proton 0.794 

Neutron 0.866 

Deuteron 0.115 

Triton 0.025 

Helion ( 3 He) 0.036 

Alpha 0.389 
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Bertini calculations for protons and neutron-induced reactions in and ^ A1 for 
1 GeV energy are given in reference 11. A direct application of equation (18) for the 
reaction yields the results given in table I. Note that these results are exact in the 
sense that the Bertini calculations are an implementation of the Serber model. The 
results in table II assume that ^Al is well approximated by an isosinglet state. Results 
for the isosinglet approximation as applied to the intermediate excited nuclear state are 
shown in table III. Based on these model calculations, the isosinglet approximation 
appears to be adequate for the description of nonelastic processes. 

CONCLUDING REMARKS 

It is interesting to see that simple algebraic relations can be derived for systems 
which react with the complexity of nonelastic reactions. In particular, when simple sys- 
tems interact with heavy target nuclei, the final reaction products do not depend strongly 
on the target isospin structure but do depend strongly on the projectile. This dependence 
appears to be in part due to the tendency of the strong force to be most attractive for sys- 
tems with small isospin numbers. The exceptional cohesiveness of states of small iso- 
spin is a property which can be of great consequence in the study of heavy -heavy nonelas- 
tic reactions. 

It is to be emphasized that these relations are not dependent on the detailed dynamics 
except to the extent that the strong force is SU(2) invariant and the isosinglet approxima- 
tion is valid. In the case of the a-^Al reaction model of Gabriel, Santoro, and Alsmiller, 
in which the model a is totally unrealistic, results that are SU(2) consistent occur. It 
should be noted that it is only the dynamical description of the a model (four unbound 
nucleons) which is unrealistic. 

From a nuclear -physics point of view, the present results yield a dynamic indepen- 
dent probe for the nonelastic process. Disagreement among secondary-particle param- 
eters for energies above that to be associated with nuclear temperature would indicate 
either SU(2) breaking terms in the first-step Serber model or more preferably violation 
of the isosinglet approximation for the residual excited nuclear state. 

The immediate consequence of the present results is to help in obtaining data for 
transport calculations. These results will reduce the number of necessary experiments 
by a factor of two. The derived relations are especially important for neutrons which are 
difficult to measure reliably and are rarely available as intense high-energy beams for 
nuclear -physics experiments. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 12, 1972. 
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because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on , the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



